Simulation of colloidal chain movements under a magnetic field 



Mareike Griinzel 

Institut fiir Theoretische Physik, Universitat zu Koln, Ziilpicher Strafle 77, D-50937 Koln, Germany 

(Dated: February 2, 2008) 

Short colloidal chains are simulated by the slithering-snake-algorithm on a simple cubic lattice. 
The dipole character of the colloidal particles leads to a long range dipole-dipole interaction. The 
solvent is simulated by the nearest neighbor Ising model. The aligning of the dipoles under a 
magnetic field gives rise to the chains to align on their part with the field direction. 
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I. INTRODUCTION 

The initial idea of this paper is based on the article of 
Goubault et al. P|, whose research partly dealt with the 
aligning of colloidal chains and their curvature behavior 
being influenced by an external magnetic field. 
The basis of the program is another program that sim- 
ulates short polymer chains with a hydrophilic head, a 
hydrophobic tail and a neutral center surrounded by oil 
and water molecules, which are simulated with the near- 
est neighbor spin Ising model • The movements of the 
chains are simulated by the slithering-snake algorithm. 
This program was created on the basis of Larson's PhD 
thesis 0. 

II. ALIGNING OF THE CHAINS 

The Ising model and the slithering-snake algorithm are 
taken over to the new program. The first aim of the 
work is to add the long-range dipole-dipole interaction 
according to the colloidal character of the chains to the 
program. The magnetic dipole interaction energy for two 
interacting magnetic momenta jli and flj (e.g. Q) is 
transformed to scalar spin variables Si = ±1 so that this 
leads to the Hamiltonian 

7i = -i ^ Juj Si Sj with 

J = / ~ • (l - 3 • (ryry )2) i, j neighbors, 

\-(.9/n,r-{l-3-(r!^/n,r) else, 

where g — jo? is the dipole strength, rf^ is the com- 
ponent of the vector rij — {rfj , rfj , rf^ ) between two re- 
garded dipoles and is the distance between them. The 
space coordinate i varies from 1 to for L x L x L lat- 
tices. The calculation of the distance is being done by 
the calculation of the coordinates of the two regarded 
dipoles with the help of the modulo-function. The x- 
coordinates, for example, are the modulo-function of the 
positions and the system width L. The difference be- 
tween the X-, y- and z-coordinates gives the distance vec- 
tor fij of the dipoles. 



Initially, the chains are aligned in the x-direction in 
every simulation. The magnetic field is oriented in z- 
direction. The computation of the energy is summed up 
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FIG. 1: Average of chain aligning of 5 chains of length 6 (con- 
centration of 0.03) in a 10 X 10 X 10 lattice at the temperature 
of 3.0 J/fcfl (100 runs). (+): x-aligning; (O): y-aligning; (*): 
z-aligning. 

for one chain interacting with the entire system due to 
the long range of the interaction. This is repeated for 
every chain. The spins of the chain molecules range from 
value Si — —1 for the tail, over to Si = for the neutral 
center, and Si — I for the second molecule over to Si = 2 
for the head. 

In this case, the magnetic field is simulated indirectly by 
aligning the dipoles in the direction of the field. The 
chains align in the field direction until equilibrium is 
achieved after approximatly 30000 time steps to be seen 
in Fig. m where the average orientations of 5 chains of 
length 6 in a 10 X 10 X 10 lattice is shown. 

The first method to calculate the aligning is the deter- 
mination of the distances of the head and the tail of the 
chains. This method favors finite-size effects in smaller 
systems, because the chains can spread over the lattice. 
- That is why only the second method has been used. It 
calculates the difference between two chain links situated 
next to each other. Then it is being decided in which di- 
rection they are standing to each other. If the difference 
is ±1, the x-aligning is raised by 1; if the difference is ±L, 
the y-aligning is raised and if it is ±L^ it concerns to the 
z-aligning. The average is calculated for every chain and 
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every run of the program at every timestep to get the 
course of chain aUgning in dependence on the time. Tak- 
ing the example of a chain with the length of 6, it can be 
arranged at a maximum of 5 in one direction. 
Simulations with a variation of the concentration of the 
chains result in a faster rise of the average alignments 
towards the field direction for small concentrations. 
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FIG. 2: Average of dipole aligning of 5 chains of length 6 
(concentration 0.03) in a 10 x 10 x 10 lattice at the temperature 
of 3.0 J/fcs (100 runs for each magnetic field strength). (+): 
x-component; (O): y-component; (*): z-component. 



III. MODIFIED MODEL 

The program has a modified version, in which the 
dipoles are three-dimensional and rotate randomly. The 
rotation takes place with the method of Barker & Watts 
in their simulation of water (1969) @. Molecules ro- 
tate with a random angle between zero and a specified 
maximum angle (in this case f) to all sides about one 
space-fixed axes chosen at random. The Hamiltonian is 
now changed as follows: 

H ~ ^ ■JiJ ^ ^ ^ ' ■^ith 

j_ - ia/rfj ) (si Sj - 3 (s, • ) {sj ■ n-i ) /rfj) 

''^ \ - [g/rlj ) {si Sj -3{si- fy )(sj ■ fy ) / r2- ) 

After the calculation of the energy to align the chains 
there is another calculation for every dipole of the chains 
relative of the rotation of the dipoles. The dipole strength 
has the value 2 for the head of every chain and 1 for the 
rest in all modified simulations. The magnetic field now 
infiuences directly the rotation of the dipoles and there- 
fore indirectly the aligning of the chains in the direction 
of the field. The magnetic field is negative in all simu- 
lations due to the fact, that a negative field affects the 



dipoles to rotate in the positive z-direction, because they 
are effective against the magnetic field. Simulations with 
a positive magnetic field result in a suppression of the 
aligning in the direction of the field. Instead, simulations 
with a negative magnetic field show the expected behav- 
ior of the dipoles under a magnetic field with various field 
strengths (Fig. [SJ. In Fig. |2lthe average value of each 
component is computed for 70 timesteps from 30 to 100 
for every magnetic field strength. This shows the equal 
arrangement of the dipoles without a magnetic field in a 
nice way, because every dipole component has the equal 
value 0.5. Simulations result in a constant value after 
30 timesteps creating a sufficient average value. The er- 
rors of these average values are too small to be visible. 
Simulations with the dipole aligning not random but in 
x-direction for the initial condition result in the same 
behavior. 
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FIG. 3: Average of the chain aligning of 5 chains of length 6 
(concentration 0.03) in a 10 x 10 x 10 lattice at the temperature 
of 3.0 J/fcs (100 runs for each magnetic field strength). (+): 
x-aligning; (O): y-aligning; (*): z-aligning. 

Fig. shows the chain aligning of 5 chains depend- 
ing on the magnetic field strength. The average value is 
computed for 850 timesteps from 150 to 1000. As Fig. 
shows, this duration is not enough to complete the align- 
ing, but its tendency is shown in Fig. For the interval 
from B = to B — 3 the aligning in the direction of 
the field does not take place due to the high tempera- 
ture of T = 3.0 J/kB- For high magnetic field strengths 
there is a saturation of the magnetic field, so that a rise 
of the aligning is no longer possible. Instead, the mid- 
dle interval of Fig. El shows a linear process of aligning 
in direction of the field. The gradient of the aligning in 
X- and y-direction is nearly the same. The fine point- 
ing into the x-direction, however, is situated a bit higher 
than the one pointing into the y-direction, because the 
chains have not arranged yet due to the short time; in 
the beginning they were arranged into the x-direction. 
Despite the small time for the average value, this is a 
clear tendency. 
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Simulations with a discrete (random angle of 7r/2 or 
— 7r/2) aligning result in the same behavior but not as 
exact as it is in the continuous aligning. The discrete 
model shows various weak points in contrast to the con- 
tinuous case. Among other things, the dipole aligning 
takes place much faster to the value 1. In relation to the 
chain aligning, there is no visible Hnear dependence on 
the magnetic field. 



pendence on the temperature or the chain length and con- 
centration can be examined. The results from Goubault 
et. al. 01 can Hkewise be researched with little effort, 
above all the curvature of the chains. The easiest way 
to examine this is taking the first method, which has 
been mentioned in section ^ For this, bigger systems 
are necessary. The Hnear dependence of the aligning on 
the magnetic field should be verified for bigger systems. 



IV. CONCLUSION 



The result is a good working program offering lots of 
possibilities to work with. Among other things, the de- 
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